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The goal of spin foam models is to provide a viable path integral formulation of 
quantum gravity. Because of background independence, their underlying framework 
has certain novel features that are not shared by path integral formulations of familiar 
field theories in Minkowski space. As a simple viability test, these features were 
recently examined through the lens of loop quantum cosmology (LQC). Results of 
\ that analysis, reported in a brief communication [1], turned out to provide concrete 

arguments in support of the spin foam paradigm. We now present detailed proofs 
of those results. Since the quantum theory of LQC models is well understood, this 
analysis also serves to shed new light on some long standing issues in the spin foam 
and group field theory literature. In particular, it suggests an intriguing possibility 
for addressing the question of why the cosmological constant is positive and small. 
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> ! I. INTRODUCTION 



Four different avenues to quantum gravity have been used to arrive at spin-foam models 
(SFMs). The fact that ideas from seemingly unrelated directions converge to the same type 
of structures and models has provided a strong impetus to the spin foam program over the 
£5 years [2]. 

The first avenue is the Hamiltonian approach to loop quantum gravity (LQG) [3-5]. 
By mimicking the procedure that led Feynman [6] to a sum over histories formulation of 
quantum mechanics, Rovelli and Reisenberger [7] proposed a space-time formulation of LQG. 
This work launched the spin-foam program. The second route stems from the fact that the 
starting point in canonical LQG is a rewriting of classical general relativity that emphasizes 
connections over metrics [8]. Therefore in the passage to quantum theory it is natural 
to begin with the path integral formulation of appropriate gauge theories. A particularly 
natural candidate is the topological B-F theory [9] because in 3 space-time dimensions it is 
equivalent to Einstein gravity, and in higher dimensions general relativity can be regarded 
as a constrained BF theory [10]. The well-controlled path integral formulation of the BF 
theory provided the second avenue and led to the SFM of Barret and Crane [11]. The third 
route comes from the Ponzano-Regge model of 3-dimensional gravity [12] that inspired Regge 
calculus in higher dimensions [13-15]. Here one begins with a simplicial decomposition of 
the space-time manifold, describes its discrete Riemannian geometry using edge lengths and 
deficit angles and constructs a path integral in terms of them. If one uses holonomies and 
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discrete areas of loop quantum gravity in place of edge lengths, one is again led to a spin 
foam. These three routes are inspired by various aspects of general relativity. The fourth 
avenue starts from approaches to quantum gravity in which gravity is to emerge from a more 
fundamental theory based on abstract structures that, to begin with, have nothing to do 
with space-time geometry. Examples are matrix models for 2-dimensional gravity and their 
extension to 3-dimensions — the Boulatov model [16] — where the basic object is a field on 
a group manifold rather than a matrix. The Boulatov model was further generalized to a 
group field theory (GFT) tailored to 4-dimensional gravity [4, 17, 18]. The perturbative 
expansion of this GFT turned out be very closely related to the vertex expansions in SFMs. 
Thus the SFMs lie at a junction where four apparently distinct paths to quantum gravity 
meet. Through contributions of many researchers it has now become an active research area 
(see, e.g., [4, 10, 19]). 

Let us begin with the first path and examine SFMs from the perspective of LQG. Recall 
that spin network states are used in LQG to construct a convenient orthonormal basis in 
the kinematical Hilbert space. A key challenge is to extract physical states from them by 
imposing constraints. Formally this can be accomplished by the group averaging procedure 
which also provides the physical inner product between the resulting states [20, 21]. From 
the LQG perspective, the primary goal of SFMs is to construct a path integral that leads 
to this physical Hilbert space. 

Heuristically, the main idea behind this construction can be summarized as follows. Con- 
sider a 4-manifold M bounded by two 3-surfaces, Si and S 2 , and a simplicial decomposition 
thereof. One can think of Si as an 'initial' surface and S 2 as a 'final' surface. One can fix a 
spin network on each of these surfaces to specify an 'initial' and a 'final' state of the quantum 
3-geometry. A quantum 4-geometry interpolating between the two can be constructed by 
considering the dual triangulation of M and coloring its surfaces with half integers j and 
edges with suitable intertwiners. The idea is to obtain the physical inner product between 
the two states by summing first over all the colorings for a given triangulation, and then over 
triangulations keeping the boundary states fixed. The second sum is often referred to as the 
vertex expansion because the M-th term in the series corresponds to a dual triangulation 
with M vertices. Since each triangulation with a coloring specifies a quantum geometry, 
the sum is regarded as a path integral over physically appropriate 4-geometries. In ordi- 
nary quantum mechanics and Minkowskian field theories where we have a fixed background 
geometry, such a path integral provides the (dynamically determined) transition amplitude 
for the first state, specified at initial time, to evolve to the second state at the final time. In 
the background independent context of quantum gravity, one does not have access to a time 
variable and dynamics is encoded in constraints. Therefore the notion of a transition in a 
pre-specified time interval is not a priori meaningful. Rather, the sum over histories now 
provides the physical inner product between solutions to the quantum constraints, extracted 
from the two spin network states. 

Over the last two years there have been significant advances in SFMs. While the structure 
of the path integral is well-motivated by the interplay between general relativity and the 
BF theory, its precise definition requires a key new ingredient — the vertex amplitude. The 
first proposal for the vertex amplitude was made over ten years ago [11]. But it turned 
out to have important limitations [22, 23]. New proposals have now been put forward 
[24-27] and, for the physically interesting regime of the Barbero-Immirzi parameter, they 
agree. Furthermore, one can regard these SFMs as providing an independent derivation of 
the kinematics underlying LQG. The detailed agreement between LQG and the new SFMs 



3 



[28, 29] is a striking development. There are also a number of results indicating that one 
does recover general relativity in the appropriate limit [32, 33]. Finally, the vertex amplitude 
is severely constrained by several general requirements which the new proposals meet. 

However, so far, the vertex amplitude has not been systematically derived following proce- 
dures used in well-understood field theories, or, starting from a well-understood Hamiltonian 
dynamics. Therefore, although the convergence of ideas from several different directions is 
impressive, a number of issues still remain. In particular, the convergence is not quite as 
seamless as one would like; some rough edges still remain because of unresolved tensions. 

For example, the final vertex expansion is a discrete sum, in which each term is itself 
a sum over colorings for a fixed triangulation. A priori it is somewhat surprising that the 
final answer can be written as a discrete sum. Would one not have to take some sort of 
a continuum limit at the end? One does this in the standard Regge approach [30] which, 
as we indicated above, is closely related to SFMs. Another route to SFMs emphasizes and 
exploits the close resemblance to gauge theories. In non-topological gauge theories one also 
has to take a continuum limit. Why not in SFMs? Is there perhaps a fundamental difference 
because, while the standard path integral treatment of gauge theories is rooted in the smooth 
Minkowskian geometry, SFMs must face the Planck scale discreteness squarely? 

A second potential tension stems from the fact that the construction of the physical inner 
product mimics that of the transition amplitude in Minkowskian quantum field theories. 
As noted above, in a background independent theory, there is no a priori notion of time 
evolution and dynamics is encoded in constraints. However, sometimes it is possible to 'de- 
parameterize' the theory and solve the Hamiltonian constraint by introducing an emergent or 
relational time a la Leibnitz. What would then be the interpretation of the spin-foam path 
integral? Would it yield both the physical inner product and the transition amplitude? 
Or, is there another irreconcilable difference from the framework used Minkowskian field 
theories? 

There is a also a tension between SFMs and GFTs. Although fields in GFTs live on 
an abstract manifold constructed from a Lie group, as in familiar field theories the action 
has a free part and an interaction term. The interaction term has a coupling constant, A, 
as coefficient. One can therefore carry out a Feynman expansion and express the partition 
function, propagators, etc, as a perturbation series in A. If one sets A = 1, the resulting 
series can be identified with the vertex expansion of SFMs. But if one adopts the viewpoint 
that the GFT is fundamental and regards gravity as an emergent phenomenon, one is led to 
allow A to run under the renormalization group flow. What then is the meaning of setting 
A = 1? Or, do other values of A have a role in SFMs that has simply remained unnoticed 
thus far? Alternatively, one can put the burden on GFTs. They appear to be efficient and 
useful calculational schemes. But if they are to have a direct physical significance on their 
own, what then would the gravitational meaning of A be? 

Such questions are conceptually and technically difficult. However, they are important 
precisely because SFMs appear to lie at a junction of several cross-roads and the recent 
advances bring out their great potential. Loop quantum cosmology (LQC) provides a physi- 
cally interesting yet technically simple context to explore such issues. In LQC the principles 
of LQG are applied to simple cosmological models which have a high degree of symmetry. 
Thanks to this symmetry, it has been possible to construct and analyze in detail quantum 
theories in a number of cases [34-47]. Furthermore, LQC shares many of the conceptual 
problems of LQG and SFMs. Therefore it provides a fertile ground to test various ideas 
and conjectures in the full theory. In the Hamiltonian context, LQC has served this role 
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successfully (for a recent review, see [48]). The goal of this paper is to first cast LQC in the 
spin foam paradigm and then use the results to shed light on the paradigm itself. 

In LQC one can arrive at a sum over histories starting from a fully controlled Hamiltonian 
theory. We will find that this sum bears out the ideas and conjectures that drive the 
spin foam paradigm. Specifically, we will show that: i) the physical inner product in the 
timeless framework equals the transition amplitude in the theory that is deparameterized 
using relational time; ii) this quantity admits a vertex expansion a la SFMs in which the 
M-th term refers just to M volume transitions, without any reference to the time at which 
the transition takes place; iii) the exact physical inner product is obtained by summing over 
just the discrete geometries; no 'continuum limit' is involved; and, iv) the vertex expansion 
can be interpreted as a perturbative expansion in the spirit of GFT, where, moreover, the 
GFT coupling constant A is closely related to the cosmological constant A. These results 
were reported in the brief communication [1]. Here we provide the detailed arguments and 
proofs. Because the Hilbert space theory is fully under control in this example, we will 
be able to avoid formal manipulations and pin-point the one technical assumption that is 
necessary to obtain the desired vertex expansion: one can interchange the group averaging 
integral and a convergent but infinite sum defining the gravitational contribution to the 
vertex expansion(see discussion at the end of section III A). In addition, this analysis will 
shed light on some long standing issues in SFMs such as the role of orientation in the spin 
foam histories [49], the somewhat puzzling fact that spin foam amplitudes are real rather 
than complex [31], and the emergence of the cosine cosS'eh of the Einstein action — rather 
than e i5EH — in the classical limit [32, 33]. 

The paper is organized as follows. In section II we summarize the salient features of 
LQC that are needed to arrive at a sum over histories formulation. Section III establishes 
the main results in the timeless framework, generally used in SFMs. In particular, we show 
that the physical inner product can be expressed as a vertex expansion. In section IV we 
introduce a deparametrization using the relational time of LQC and obtain an equivalent but 
distinct vertex expansion, more directly related to the transition amplitude. The existence 
of distinct vertex expansions which sum to the same result suggests the possibility that there 
may well be distinct but physically equivalent vertex amplitudes in SFMs, each leading to a 
perturbative expansion that is tailored to a specific aspect of the physical problem. To avoid 
repetition, we adopted a strategy that is opposite of that used in [1]: here we provide detailed 
derivations in the timeless framework (section III) and leave out the details while discussing 
analogous results in the deparameterized picture (section IV). Section V summarizes the 
main results and discusses some generalizations and open issues. A number of technical 
issues are discussed in three Appendices. 

II. LQC: A BRIEF OVERVIEW 

We will focus on the simplest LQC model that has been analyzed in detail [34-36, 39]: 
the k=0, A=0 Friedmann model with a massless scalar field as a source. However, it should 
not be difficult to extend this analysis to allow for a non-zero cosmological constant [40, 41] 
or anisotropies [43, 44] or to the spatially compact k=l case [37]. 

In the FRW models, one begins by fixing a (spatial) manifold S, topologically M 3 , Carte- 
sian coordinates x % thereon, and a fiducial metric q° b given by q° b dx a dx b = dx\ + dx\ + dx\. 
The physical 3-metric q a b is then determined by a scale factor a; q a b = a 2 q° b . For the 
Hamiltonian analysis one fixes a cubical fiducial cell V whose volume with respect to q° ah is 
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V Q so that its physical volume is V = a 3 V a . The quantity v defined by V = 2n / -f£p 1 \v\ turns 
out to be a convenient configuration variable, where 7 is the Barbero-Immirzi parameter of 
LQG [39]. 1 

The kinematical Hilbert space is a tensor product "H kin = "H^^^kin" °f the gravitational 
and matter Hilbert spaces. Elements ^{y) of "Hf[„ v are functions of v with support on a 
countable number of points and with finite norm 

ll^ll 2 : = Z)„ \^( v )\ 2 - The matter Hilbert 
space is the standard one: "Hj^f* = L 2 (R, d(f>). 2 Thus, the kinematic quantum states of the 
model are functions \l/(z/, 0) with finite norm ||^/|| 2 := J2 U Jd0 [^(z/, 0)| 2 . A (generalized) 
orthonormal basis in "H kin is given by \u, 0) with 

(^0>,0) = <W<5(0'>)- (2.1) 

To obtain the physical Hilbert space, one first notes that the quantum constraint can be 
written as 

- C*(v, 0) = 0) + 0) = (2.2) 

where 6 is a positive and self-adjoint operator on [50]- More explicitly, 6 is a second 

order difference operator [43] 

(6* ) (1/) := [ y/\v{y + U )\ {v + 2Q *(v + AQ - 2v 2 ^{y) 



+ - U )\{v -2Q^{v - U )] , (2.3) 



where £ Q is related to the 'area gap' A = 4v / 3tt7^pi via t 2 Q = A. The form of shows that 
the space of solutions to the quantum constraint can be naturally decomposed into sectors 
in which the wave functions have support on specific '//-lattices' [35]. For definiteness, we 
will restrict ourselves to the lattice v = 4n£ where n is an integer. Details of the expression 
of will not be needed in most of our analysis. 

The scalar field is monotonic on all classical solutions (also in the cases when k=l, 
and A 7^0) and therefore serves as a relational time variable, a la Leibnitz, in the classical 
theory. This interpretation carries over to the quantum theory. For, the form of the quantum 
constraint (2.2) is similar to that of the Klein-Gordon equation, playing the role of time 
and —0 of the spatial Laplacian (or, the elliptic operator generalizing the Laplacian if we are 
in a general static space-time). Therefore, in LQC, one can use as an internal time variable 
with respect to which physical quantities such as the density, scalar curvature, anisotropies 
and shears in the Bianchi models [43, 44], and the infinitely many modes of gravitational 
waves in the Gowdy models [45-47], evolve. 

In the spin foam literature, by contrast, one does not have access to such a preferred 
time and therefore one chooses to work with the timeless formalism. Therefore let us first 



1 In LQG the basic geometric variable is an orthonormal triad and the physical metric q a b is constructed from 
it. If the triad has the same orientation as the fiducial one, given by the coordinates x l , the configuration 
variable v is positive and if the orientations are opposite, v is negative. Physics of the model is insensitive 
to the triad orientation and hence to the sign of v. In particular the kinematic and physical quantum 
states satisfy SS?(v,<j)) = ^>{—v,4>). 

2 One can also use a 'polymer quantization' of the scalar field at the kinematical level but the final physical 
theory turns out to be the same. 
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forgo the emphasis on using as internal time and simply implement the group averaging 
procedure which uses the constraint operator as a whole, without having to single out a 
preferred time variable [20, 21]. This procedure plays an important role in sections III and 
IV. Therefore it is useful to summarize it in some detail. One begins by fixing a dense 
sub-space S of "Hkin- In LQC, this is generally taken to be the Schwartz space of smooth 
functions f(u, 0) which fall off to zero at infinity faster than any polynomial. The first step 
in the group averaging procedure is to extract a solution ^ff(u, 0) to the quantum constraint 
operator (2.2) from each / G S. These solutions are not normalizable in "H kin because the 
spectrum of the constraint C on "H kin is continuous. The second step of the group averaging 
procedure provides an appropriate inner product between solutions ^f(v, <f>)- 

Denote by e k (u), with k G (—00,00) a complete set of orthonormal eigenfunctions of 
on H^iii ' ■ We will denote the eigenvalues by io\ and, without loss of generality, assume that 
u k > [35, 36]. (Eigenfunctions and operator functions of are discussed in Appendix C.) 
Any f(u, 0) G S can be expanded as 

f(v,4>) = fdk ±jd P4> KKj>t)j*+*e k (y). (2.4) 

Here and in what follows the range of integrals will be from —00 to 00 unless otherwise 
stated. Using this expansion, we can group-average any f{y, 0) to obtain a distributional 
solution (in S*) \I//(z/, 0) to the quantum constraint: 

¥/(!/, := Jda[ e - c 2b |/(z/,0)] = Jdfc /dp, <S(p} - a# 2\ P4> \f{k^) e*** e k (u) , (2.5) 

where, the operator 2|p^| has been introduced just for later technical simplification. Had 
we dropped it, we would have associated with / the solution (2|p^|) _1 tyf and, in the end, 
obtained a unitarily equivalent representation of the algebra of Dirac observables. 

By carrying out the integral over the expression of \&/ can be brought to the desired 
form: 

*/(!/, 0) = Jdk[f(k,u k )e^e k (v) + f(k,-u k )e-^e k (u)} 

=: *+(!/, 0) + 0). (2.6) 

By their very definition ^^(1/, 0) satisfy 

^(z/,0) = e ± ^^°)^(z/,0 o ), (2.7) 

whence they can be interpreted as 'positive and negative frequency solutions' to (2.2) with 
respect to the relational time 0. Thus the group average of / is a solution tyf to the quantum 
constraint (2.2) which, furthermore, is naturally decomposed into positive and negative 
frequency parts, "iff is to be regarded as a distribution in S* which acts on elements g G S 
via the kinematic inner product [20, 21]: 

= fdk Jdp <j> 5(p 2 4) -ujl)2uj k f(k,p < f > )g(k,p < f > ) 

= Jdk [f(k,u k )g(k,u k ) + f(k, -u k )g(k, -w k )] . (2.8) 
Finally, the group averaged scalar product on solutions is given just by this action 



7 



[20, 21]. Thus, given any elements /, g in S, the scalar product between the corresponding 
group averaged states ^ f,^ g is given by 

(¥/, * g ):=(? f \g)=Wg\f). (2.9) 

In section III we will obtain a vertex expansion for this scalar product. 

A conceptually important observation is that, as in the Klein-Gordon case, there is a 
superselection. A complete set of Dirac observables is given by the scalar field momentum 

= —id^ and the volume V\^ a (or, equivalently, the energy density operator p\ ( f )o ) at the 
value = O of the internal time. (The factor of \p^\ introduced above simplifies the explicit 
expressions of V\ ( p and p\^, [35, 36, 39].) The action of these Dirac observables as well 
as time evolution leaves the space of positive and negative frequency solutions invariant. 
Therefore, as in the Klein-Gordon theory, we are led to work with either set. In LQC, one 
generally works with the positive frequency ones. Then the physical Hilbert space "H p h y of 
LQC consists of positive frequency solutions \l/ + (z/, 0) to the quantum constraint (2.2), i.e. 
solutions satisfying 

-id+ 0) = VQV+{v, 0) = if 0) (2.10) 
with inner-product (2.9). This inner product can be re-expressed simply as: 

(* + , $ + ) phy = £ *+(!/, (2.11) 

u=4n£ 

and is independent of the value O of at which the right side is evaluated. 

While this construction of "H p h y does not require us to think of as internal time in 
quantum theory, this interpretation is natural in the light of final Eqs (2.10) and (2.11). 
For, these equations suggest that we can think of v as the sole configuration variable and 
introduce 'Schrodinger states' ^{y) through the physical inner product (2.11). These 'evolve' 
via (2.10). This is the 'deparameterized' description to which we will return in section IV. In 
this picture, the restriction to positive frequency states has direct interpretation: p^ = \/0 
is now a positive operator on % p h y just as p is a positive operator on the traditional Klein- 
Gordon Hilbert space. 



III. THE TIMELESS FRAMEWORK 

Recall that in the spin foam literature, one works with the timeless framework because 
a natural deparametrization is not available in general. To mimic the general spin foam 
constructions in LQC, in this section we will largely disregard the fact that the scalar field 
can be used as relational time and that the final constraint has the form of the Schrodinger 
equation. Instead, we will use the group averaging procedure for the full constraint 

c = -di-e = P l-e (3.i) 

and incorporate the positive frequency condition in a second step. None of the steps in this 
analysis refer the evolution in relational time mentioned above. Thus, the primary object 
of interest will be the physical scalar product, rather than the transition amplitude for a 
Schrodinger state \I/(z/, 0;) at an initial 'time instant' 0« to evolve to another state $(z/, 0/) 
at a final 'time instant' 0/. 
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In section II we considered general kinematic states f{v,fa). In this section, by contrast, 
we will focus on the basis vectors \u, fa) in "Hkm which are the LQC analogs of spin networks 
that are used to specify the boundary states in SFMs. Following the setup introduced in 
section I let us then fix two kinematic states, \vi,fa) and \i/f,(f>f). For notational simplicity, 
we will denote the group averaged solutions to (2.2) they define by |[z/j,0j) and |[^/, </>/])• 
The group averaged inner product between these states is given by 

{[vfAfl [viAi\) = lS&a(v h (t>f\e iaC MK0*>. (3.2) 

Our goal is to express this scalar product as a vertex expansion a la SFMs and study 
its properties. In section III A we will begin by rewriting it as a sum over histories a la 
Feynman [6] and then rearrange the sum as a vertex expansion. In section III B we will 
arrive at the same expansion using perturbation theory in a suitably defined interaction 
picture. This procedure is reminiscent of the perturbation expansion used in GFTs. As 
an important consistency check, in section III C we verify that this perturbative expansion 
does satisfy the constraint order by order. Finally, in section HID we observe that, in 
this simple example, the coupling constant A used in the expansion is intimately related to 
the cosmological constant A. Although the precise relation we obtain is tied to LQC, the 
observation illustrates in a concrete fashion how one may be able to provide a gravitational 
interpretation to A in GFTs and suggests an avenue for GFT to account for the smallness 
of A. 

A. Sum over Histories 

Following Reisenberger and Rovelli [7], let us first focus on the amplitude 

A(i/f, 0/; v h fa; a) = 2 (is f , <f> f \ e taC \p^\ \v h fa) (3.3) 

which constitutes the integrand of (3.2). Mathematically one can choose to regard aC 
as a Hamiltonian operator. Then A(vf, 0/, v { , fa, a) can be interpreted as the probability 
amplitude for an initial kinematic state \vi,fa) to evolve to a final kinematic state \vf,(f>f) 
in a unit 'time interval' and we can follow Feynman's procedure [6] to express it as a sum 
over histories. Technically, a key simplification comes from the fact that the constraint C is 
a sum of two commuting pieces that act separately on "Hj^f 1 and "Hf[„ v - Consequently, the 
amplitude (3.3) factorizes as 

A(v f , 0/; Ui, fa; a) = A^fa, fa; a) A G {v s , v { ; a) (3.4) 

with 

Mfa, fa; a) = 2 (<P f \e iap l \ P <p\\fa), and A G (is f , u t ; a) = (u f \e' tae \u t ) . (3.5) 

It is easy to cast the first amplitude, A^, in the desired form using either a standard Feynman 
expansion or simply evaluating it by inserting a complete eigen-basis of p^. The result is: 

A^f, fa; a) = 2 / d P(j> e l < e **(*/-&) \ v< ^\ (3.6) 

The expansion of the gravitational amplitude Aq is not as simple. We will first express 
it as a sum over histories. In a second step, we will evaluate the total amplitude (3.3) by 
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integrating over a for each history separately. Although it is not a priori obvious, we will find 
that the amplitude associated to each history is manifestly finite and the total amplitude 
can be written as a discrete sum that mimics the vertex expansion in SFMs. 

1. The gravitational amplitude Aq 

As mentioned above, to apply the standard Feynman procedure we will regard e~ ta@ as 
an 'evolution operator' with 'Hamiltonian' a© and a 'time interval' At = 1. We emphasize 
that this 'evolution' is a just a convenient mathematical construct and does not correspond 
to the physical evolution with respect to the relational time variables normally used 
in LQC. Rather, since it is generated by the constraint C, physically it represents gauge 
transformations (or time reparameterizations). 

Let us divide the interval At = 1 into N parts each of length e = 1/N and write the 
gravitational amplitude Aa(vf, vf, a) as 

(v f \e- ta@ \v t ) = J2 (^|e" l£Q0 |^-i)(^-i|e-" ae |^_ 2 )...(^|e- l£Qe |z/,) (3.7) 

PjV-l,...,Pl 

where we have first split the exponential into N identical terms and then introduced a 
decomposition of the identity operator at each intermediate 'time' r = ne, n = 1,2, .., N — 1. 
For notational simplicity, we will denote the matrix element {v n \e~ teae \v n -i) by Up n p n _ l and 
set Uf = u N and z/j = u Q . We then have 

A^/, v%\ ex) — ^2 ^atPjv-i U Vn _ iVn _ 2 . . . U PlPo . (3.8) 

PjV_i,...,Pl 

The division of At provides a skeletonization of this 'time interval'. An assignment a N = 
(pjv, • • • , vq) of volumes to the N + 1 time instants t = en can be regarded as a discrete 
(gauge) history associated with this skeletonization since one can envision the universe going 
from v n -i to v n under a finite 'evolution'. The matrix element is given by a sum of amplitudes 
over these discrete histories with fixed endpoints, 

A G {v h Vi\ a) = ^ A(a N ) = U UnPn _ 1 U Dn _ i9n _ 2 . . . U V2Vl U PlPo . (3.9) 

(Tjv &N 

The next step in a standard path integral construction is to take the 'continuum' limit, 
N — > oo, of the skeletonization. In particle mechanics at this stage one uses a continuous 
basis (say the position basis \x)) to carry out this expansion. By contrast, our basis \u n ) 
is discrete. As a result, one can make rigorous sense of the iV — > oo limit by reorganizing 
the well-defined sum (3.9) according to the number of volume transitions. The remainder 
of section III A 1 is devoted to carrying out this step. 

This task involves two key ideas. Let us first note that along a path o^r, the volume v 
is allowed to remain constant along a number of time steps, then jump to another value, 
where it could again remain constant for a certain number of time steps, and so on. The 
first key idea is to group paths according to the number of volume transitions rather than 
time steps. Let us then consider a path which involves M volume transitions (clearly, 
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M < N): 

N 2 

a N — ( v Mi ■ ■ ■ , vm] vm-\, ■ ■ ■ , vm-\\ ; V\, ■ ■ ■ i V\\ Ifo, ■ ■ • , vq) ■ (3.10) 

Thus, the volume changes from v m -i to u m at 'time' r = N m e and remains v m till time 
t = N m+ i e. Note that v m is distinct from D m used in (3.9): While v m is the volume after 
the m-th volume transition along the given discrete path, v m is the volume at the end of 
the m-th time interval, i.e., at r = me. 

These discrete histories can be labeled more transparently by two ordered sequences 

a N = { ( v m, vm-i, ■ ■ ■ , V\, u ); (N M ,N M -i, ■ ■ ■ ,N 2 ,Ni)}, v m ^ v m -i, 

(3.11) 

where v M , . . . , v denote the volumes that feature in the history cxjjf and N k denotes the 
number of time steps after which the volume changes from Uk-i to vy.. Note that while no 
two consecutive volume values can be equal, a given volume value can repeat in the sequence; 
v m can equal some v n if n ^ m ± 1. The probability amplitude for such a history is 
given by: 

^(°"tv ) = \Uv m vm\ N Nm 1 ^v m vm-i ■ ■ ■ [^lfi]^ 2 Nl 1 U V1VQ [U U0UQ ] Nl 1 . (3-12) 

The second key idea is to perform the sum over all these amplitudes in three steps. First 
we keep the ordered set of volumes (um, • • • , ^o) fixed, but allow the volume transitions to 
occur at any value r = ne in the interval Ar, subject only to the constraint that the m-th 
transition occurs before the (m+l)-th for all m. The sum of amplitudes over this group of 
histories is given by 

N-l JVm-1 n 2 -i 
A N (u M ,...,u ;a)= E ■■■J2M°n)- (3-13) 

N M =M N M -i=M-l ATi=l 

Next we sum over all possible intermediate values of v m such that v m 7^ v m -ii keeping 
vq = Vi, u M = Vf to obtain the amplitude A N (M) associated with the set of all paths in 
which there are precisely M volume transitions: 

A N (M;a)= A N (v M ,. . .,v ;a) (3.14) 

VM-1, — ,V1 

Finally the total amplitude Aci^f] Vi, a) is obtained by summing over all volume transitions 
that are permissible within our initially fixed skeletonization with iV time steps: 

N 

A G (v f , a) = J2 a n{M; a) (3.15) 

M=0 

This concludes the desired re-arrangement of the sum (3.9). The sum on the right side is 
manifestly finite. Furthermore, since Ac{vf, vu a ) — { v f\ e ~ Wi), the value of the amplitude 
(3.15) does not depend on N at all; the skeletonization was introduced just to express this 
well-defined amplitude as a sum over histories. Thus, while the range of M in the sum and 
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the amplitude A^(M; a) in (3.15) both depend on N, the sum does not. 

Therefore we are well positioned to get rid of the skeletonization altogether by taking 
the limit N goes to infinity. Note first that with our fixed skeletonization, the gravitational 
amplitude is a finite sum of terms, 

Aa(u f , Vi\ a) = A N (0; a) + A N (1; «) + ... + A N (M; «) + ... + A N (N; a) (3.16) 

each providing the contribution of all discrete paths that contain a fixed number of volume 
transitions. Let us focus on the Mth term in the sum: 

A N (M;a)= ^ A N (u M , . . . ,u ;a) (3.17) 

Now, in Appendix A we show that the limit liniAr^oo An(um, ■ ■ ■ , v ; ot) exists and is given 
by 

A(iy M ,...,u ;a) : = lim A N (is M , . . . , i/ ; a) 

N— >oo 

= Jldr M jQ M drA f _i ...J^dn A(u M , . . . , v ; t m , ■ ■ ■ , n; a) (3.18) 

where 

A(v M , . . • , z/ ; r M , • • • , ri; a) := e ~ i(:i -™ )aB w (-iaQ^^J x 

... e-i^-njae^^ (-iae^) e _<naG ^ . (3.19) 

Note that the matrix elements ® Vm v n = (^ml©!^™) of in "Hf[„ v can be calculated easily 
from (2.3) and vanish if {y m — v n ) ^ {0, ±4£ }- Therefore, explicit evaluation of the limit is 
rather straightforward. We will assume that the limit iV — > oo can be interchanged with the 
sum over vm-i, ■ ■ ■ v\- (This assumption is motivated by the fact that in the expression of 
A{i>Mi . . . ,vq;cx) most matrix elements of vanish, and since the initial and final volumes 
are fixed, the sums over intermediate volumes vm-i, ■ ■ ■ , v\ extend over only a finite number 
of non-zero terms.) Then it follows that 

A G (M;a) := lim A N (M;a) 

N->oo 

exists for each finite M. Note that the reference to the skeletonization disappears in this 
limit. Thus, Aa(M;a) is the amplitude obtained by summing over all paths that contain 
precisely M volume transitions within the 'time interval' At = 1, irrespective of precisely 
when and at what values of volume they occurred. Finally, (3.16) implies that the total 
gravitational amplitude can be written as an infinite sum: 

oo 

A G (u f , Vi ; a)=J2 MM; a) (3.20) 

M=0 

While each partial amplitude Aq(M; a) is well-defined and finite, it does not ensure that 
the infinite sum converges. A priori the infinite sum on the right hand side of (3.20) could 
be, for example, only an asymptotic series to the well-defined left side. Also, our derivation 
assumed that the limit N — > oo commutes with the partial sums. Both these limitations will 
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be overcome in section III B: We will see that A G (vj, uf, a) is indeed given by a convergent 
sum (3.20). 

The expression (3.18) still contains some integrals. These can be performed exactly. The 
case when all of (z/ M , . . . , u ) are distinct is straightforward and the result as given in [1]. 
The general case is a little more complicated and is analyzed in Appendix B. The final 
result is: 



A(vm, ■ ■ ■ , Z^O; a ) — ®v m vm-i®vm-ivm-2 ■ ■ ■ ®^2fi®flfo X 

P 1 / Q X"*- 1 P e -i a e WmWm A T 

11 ( n , _iv ( ) 1^t\p~7q ie ) ^ 3 - 21 ^ 

where, since the volumes can repeat along the discrete path, w TO label the p distinct values 
taken by the volume and n m the number of times that each value occurs in the sequence. 
The n m satisfy n\ + . . . + n p — M + 1. 

To summarize, we have written the gravitational part Ao(vf, vf, a) of the amplitude as a 
'sum over histories': 

oo 

Aa{v f , Vi\ a) = ^ ^2 A(v M , . . . , v ; a) (3.22) 

M=0 v m-i,—,v\ 

with A(u M , . . . , f ; a) given by (3.21). This expression consists of a sum over M, the number 
of volume transitions, and a sum over the (finite number of) sequences of M — 1 intermediate 
volumes that are consistent with the boundary conditions and the condition that v m ^ v m +i- 
In section III A 2 we will use this sum to generate the 'vertex expansion' of the physical inner 
product. 



2. Vertex expansion of the physical inner product 

Recall that the group-averaged scalar product can be expressed as 

([zz/,0/], [u i ,<j> i ])=2jdaA^(u i ,<l> i ;a)A G (i'f,u i ;a). (3.23) 

The main assumption in our derivation — the only one that will be required also in section 
III B — is that one can interchange the integration over a and the (convergent but infinite) 
sum over M in the expression of Ac(vf, Vi\ a). Let us then use expressions (3.6) and (3.22) 
of A^ and Aq, make the interchange and carry out the integral over a. The scalar product 
(3.23) is then re-expressed as a sum of amplitudes associated with discrete paths (vm, ■ ■ ■ ,^o)' 

oo 

vm, ■ ■ ■ ,vo] 4>f, <pi) , (3.24) 

M=0 "m-1i-i"i 
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where, 

A{U M , ■ ■ ■ , VQ\ (fif, <Pi) = 2 QumVM-1 &VM-IVM-2 ■ ■ ■ x (3.25) 
P 1 / Q \ P 

Ilr^nT ( TxV— ) E /dp*^/-W M n /y— ) • 

The right side is a sum of distributions, integrated over p^. It is straightforward to perform 
the integral and express A(i/m, ■ ■ ■ ,^o] (pf, (pi) m terms of the matrix elements of 0: 

A(U M , • • • , Z^O; 0/, 4>i) = ®VM-lVM-2 ■ ■ ■ ^Wl ©^o x (3.26) 

where A(p — (pf — (pi- Since by inspection each amplitude A(i/ M , . . . , u , <f>f, (pi) is real, the 
group averaged scalar product (3.24) is also real. 

Finally, as explained in section II, the group averaging procedure yields a solution which 
has both positive and negative frequency components while the physical Hilbert space con- 
sists only of positive frequency solutions. Let us denote the positive frequency parts of the 
group averaged ket \[is, <p]) by \[u, (p]+)- Then, the physical scalar product between these 
states in Uphy is given by a sum over amplitudes A(M), each associated with a fixed number 
of volume transitions: 

(WfAfU, h,0i] + )phy 



X 

(Note that the right side is in general complex, a point to which we will return in section 
V.) This is the vertex expansion of the physical inner product we were seeking. It has 
two key features. First, the integral over the parameter a was carried out and is not 
divergent. This is a non-trivial and important result if we are interested in computing 
the physical inner product perturbatively, i.e., order by order in the number of vertices. 
Second, the summand involves only the matrix elements of which are easy to compute. 
As remarked earlier, significant simplification arises because Eq (2.3) implies that UmUn is 
zero if v m -v n ^ {0, ±4£ }. 

Let us summarize. We did not begin by postulating that the physical inner product is 
given by a formal path integral. Rather, we started with the kinematical Hilbert space and 
the group averaging procedure and derived a vertex expansion of the physical inner product. 
Because the Hilbert space framework is fully under control, we could pin-point the one 
assumption that is needed to arrive at (3.27): the sum over vertices and the integral over a 
can be interchanged. In the full theory, one often performs formal manipulations which result 
in divergent individual terms in the series under consideration. (For instance sometimes one 



E A(M) (3.27) 

M=0 

oo 

^M"Jlf-l ®V M _ 1 V M _ 2 ■ ■ ■ ®v 2 Ui 01/1^0 

M=0 v M-\,—,v\ 



m + 1 



P I , Q y k -l p e iy/B WmWm A<f> 

11 ( n , - TTi v^ft ) 2^TW 



k=1 ( n k ^y-^-dQ WkWk j ! n^^wffl QwjWj)- 
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starts by expanding the very first amplitude (3.3) in powers of a even though the a integral of 
each term is then divergent [7, 19]). In our case, individual terms in the series are all finite, 
and, as we will show in section III B, even the full series (3.22) representing the gravitational 
amplitude is convergent. Nonetheless, at present the interchange of the a-integral and the 
infinite sum over M has not been justified. If this gap can be filled, we would have a fully 
rigorous argument that the well-defined physical inner product admits an exact, convergent 
vertex expansion (3.27). (This assumption is needed only in the timeless framework because 
the integration over a never appears in the deparameterized framework of section IV.) In 
particular, there is no need to take a 'continuum limit'. 



B. Perturbation Series 

We will now show that the expression (3.27) of the transition amplitude can also be 
obtained using a specific perturbative expansion. Structurally, this second derivation of the 
vertex expansion is reminiscent of the perturbative strategy used in group field theory (see, 
e.g., [17, 18]). 

Let us begin by considering the diagonal and off-diagonal parts D and K of the operator 
6 in the basis \v = 4n£ ). Thus, matrix elements of D and K are given by: 

D v ' v — Q uv 8 v i v , K v i v = \^ vv v '^t= v (3.28) 

Clearly C = — D — K. The idea is to think of — D as the 'main part' of C and K as 
a 'perturbation'. To implement it, introduce a 1-parameter family of operators 

C x =pl-e x :=pl-D-\K (3.29) 

as an intermediate mathematical step. The parameter A will simply serve as a marker to 
keep track of powers of K in the perturbative expansion and we will have to set A = 1 at 
the end of the calculation. 

Our starting point is again the decomposition (3.4) of the amplitude A(vf,<f>f;i/i,<f>i;a) 
into a scalar field and a gravitational part. The A dependance appears in the gravitational 
part: 

Ag\v f ,Vi,a) := (u f \e-^\u t ). (3.30) 

Let us construct a perturbative expansion of this amplitude. Again we think of e~ ta@x as 
a mathematical 'evolution operator' defined by the 'Hamiltonian' aQ\ and a 'time interval' 
At = 1. The ' unperturbed Hamiltonian' is aD and the 'perturbation' is XaK. Following 
the textbook procedure, let us define the 'interaction Hamiltonian' as 

H^r) = e iaDT aK e~ iaDl ". (3.31) 

Then the evolution in the interaction picture is dictated by the 1-parameter family of unitary 

/grav 
kin 

U x (t) = e iaDr e- iae ' T , satisfying dU ^ T > = -iXH^Ux^) . (3.32) 

dr 



operators on 
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The solution of this equation is given by a time-ordered exponential: 
U x (r) = Te- l f T o H ^ dr 

oo 

= E XM /o dr M f M dr M ^ . . . ftdn [-iH^tm)] ... HHjin)] . (3.33) 



M=0 



Next we use the relation e~ iaex = e~' iaD U\(l), with U\ given by (3.33), take the matrix 
element of e ia0A between initial and final states, = z/ ) and \uf = vm), and write out 
explicitly the product of the if/'s. The result is 

oo 

Ag\u f ,u i ,a) = '£ X M Jldr M ...^dr 1 Y, VM _ u ..., Vl b~ i{1 ~™ )aDvMVM \ * 



M=0 



-iaK VMVM _ x ) ... (-taK UlU0 ) [e^"] (3.34) 



We can now replace D and K by their definition (3.28). Because K has no diagonal matrix 
elements, only the terms with v m ^ v m +i contribute and the sum reduces precisely to 



A 



oo 

(A) 
G 



(v f , Ui, a) = E XM [ E A{yMi ...,v ;a) , (3.35) 



M=0 vm-i,—,vi 



where A(vm, ■ ■ ■ , W, ot) is given by (3.21) as in the sum over histories expansion of section 
III A 1. 

We can now construct the total amplitude by including the scalar field factor (3.6) and 
performing the a integral as in section III A 2. Then the group averaged scalar product is 
given by 

oo 

([^,0/], K<i>i]) W = Y,\ M [ J2 ^m,...,^o,0/,0*)] (3.36) 

M=0 ^M-lv-! !/ l 

where A(vm, ■ ■ ■ , vo, 0/, <f>i) is given in (3.26). If we now set A = 1, (3.36) reduces to (3.24) 
obtained independently in section III A 2. 

Finally, let us restrict ourselves to the positive frequency parts <f>] + ) of [v,(f>]) which 
provide elements of "H p h y . Reasoning of section III A 2 tells us that the physical scalar 
product ([vf,(f>f]+, [vi,(/)i\ + ) phy is given by (3.27). 

Thus, by formally regarding the volume changing, off-diagonal piece of the constraint 
as a perturbation we have obtained an independent derivation of the vertex expansion for 
([i/f, 0/] + , [vi, 0i]+) p hy as a power series expansion in A, the power of A serving as a bookmark 
that keeps track of the number of vertices in each term. In this sense this alternate derivation 
is analogous to the vertex expansion obtained using group field theory. This derivation has 
a technical advantage. Since Hi is self-adjoint on Hf™^ , it follows that the expansion (3.33) 
of U\(t) is convergent everywhere on H^m [51]- This in turn implies that the right hand 
side of (3.35) converges to the well-defined gravitational amplitude A^q = {vf\e~ ia@x \i>i) . 
However, to arrive at the final vertex expansion starting from (3.35) we followed the same 
procedure as in section III A 2. Therefore, this second derivation of the vertex amplitude 
also assumes that one can interchange the integral over a with the (convergent but) infinite 
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sum over M in (3.35). 

C. Satisfaction of the constraint 

The physical inner product between the basis states defines a 2-point function: 

G(v f ,(f) f ; vrfi) := ([*//, 0/]+, h,0i]+) p hy (3.37) 

and it follows from section II that it satisfies the constraint equation in each argument. Since 
G(uf, (fif] z/j0j) = G(vi, (pi, Vf(f>f), it suffices to focus just on one argument, say the final one. 
Then we have: 

[% f -e f ]G(v f ,(i> f ;v i ,<f> i ) = (3.38) 

where 0/ acts as in (2.3) but on Vf in place of v. If one replaces G by 0\, one obtains a 
2-point function G\(vf,<f>f; Vi,(f>i) which, as we saw in section III B admits a perturbative 
expansion: 

oo 

Gx(is f ,<j)f, Vi,pi) = ^ ^ M ^M(^/,0/;^,0i), (3.39) 

M=0 

where Am is the amplitude defined in (3.27): 

A M (v f ,(j) f ]v h (j)i) = ^2 a +(vm, ■ ■ .fo; (t>fAi) 

vm-i,---,v\ 



fM-l,-,"l 



The suffix + in A + (um, . . . , z/ ; 0j, 0j) emphasizes that we have taken the positive frequency 
part. 

As a non-trivial check on this expansion we will now show that G\ satisfies (3.38) order 
by order. Since Q\ = D + \K, our task reduces to showing 

(&l f - D f ) A M {y h 4> f ; Ui<f>i) - K f A M -i(vf, </)f, Mi) = . (3.41) 

We will show that the left hand side is zero path by path in the sense that for every path 
acted on by the off-diagonal part there are two paths acted on the diagonal part that cancel 
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it. Without loss of generality we assume that Vf = w p in (3.40). Then we have 

..0, 



( d l f - D f) A +{ V fi "M-l, ■ ■ ■ , V\, V%\ • 



n (n k 

,k=l K ft 



d 



n k -l P 



1 



d 



w k w k . 

n k -l P 



E 



e 



w m w m '- 



t A0 



'w k w k 



(3.42) 



If Wp occurs with multiplicity n p — 1, if Vf is the only volume to take the value w p then 
there are no derivatives in Q WpWp in the above equation and it simplifies to 



( 9 l f - D f) A +( U f> "M-l, Vi\ 0/, <t>i) = ®v f VM-i Q VM-lVM-2 



■ ■ ■ U2 uiO ulUi x 



p-i 



P {n k _ i)i \ae WkWk 



d 



n k -l 



E 

m=l 



rijyjl^WmWm ®WjWj) 



= OvfUM^A+il/M-l, ■ ■ ■ , V X , Vi\ 0/, 0i) . (3.43) 

Thus, on simple paths where the final volume occurs only once in the sequence, the action 
of [dj, — D] is to give the amplitude of the path without Vf, times a matrix element of 
9 related to the transition from vm-i to v$. In general, the value of the final volume can 
be repeated in the discrete path; n p ^ 1. In that case we need to push Q WpWp under the 
derivatives but the final result is the same. Thus, in all cases we have 

( d l f ~ D f) A +("f, "M-i, ■ ■ ■ , vi, v%\ 0/, <l>i) = ®vf M -i a +(vm-i, Vi\ 0/, (f>i) • (3.44) 

Finally, it is straightforward to evaluate the action of the off-diagonal part on Am-i (see 
(3.41)): 

KA+(vf, v M -2, ■ ■ ■ v%\ (pf,(pi) = ^2 ®vfM-iA+{vM-i, vm-2, ■■■,v\, v%\ (pf,(pi) • (3.45) 

v M-\ 

Combining these results we see that Eq. (3.41) is satisfied. Thus the vertex expansion we 
obtained is a solution to the quantum constraint equation. Further it is a good perturbative 
solution in the sense that, if we only take paths in which the number of volume transitions 
is less than some M*, then the constraint is satisfied to the order X M : 

M* 

[dl f - (D f + XK f )] A M A M (u f , /5 u it 0») = 0(X M ^) (3.46) 

M=0 

Also in this calculation the cancelations occur in a simple manner; the off-diagonal part 
acting on paths with M — 1 transitions gives a contribution for each path with M transitions 
that could be obtained by a adding a single additional transition in the original path. These 
contributions cancel with the action of the diagonal part on the paths with M transitions. 
This calculation provides an explicit check on our perturbative expansion of the physical 
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inner product. This is a concrete realization, in this simple example, of a central hope of 
SFMs: to show that the physical inner product between spin networks, expressed as a vertex 
expansion, does solve the Hamiltonian constraint of LQG order by order. 



D. The 'coupling constant' A and the cosmological constant A 

So far we have regarded the GFT inspired perturbation theory as a calculational tool and 
the coupling constant A as a book-keeping device which merely keeps track of the number 
of vertices in the vertex expansion. From this standpoint values of A other than A = 1 have 
no physical significance. However, if one regards GFT as fundamental and gravity as an 
emergent phenomenon, one is forced to change the viewpoint. From this new perspective, 
the coupling constant A is physical and can, for example, run under a renormalization group 
flow. The question we raised in section I is: What would then be the physical meaning of A 
from the gravitational perspective? Surprisingly, in the LQC model under consideration, A 
can be regarded as (a function of) the cosmological constant A. 

Let us begin by noting how the quantum constraint changes in presence of a cosmological 
constant A: 

-C(A) = &l + 0(A) = &l + 6 - 7rG 7 2 Az/ 2 . (3.47) 

Thus, only the diagonal part of G is modified and it just acquires an additional term propor- 
tional to A. In the GFT-like perturbation expansion, then, we are led to decompose 0a(A) 
as 

A (A) = -D(A) + XK where D(A) = ttG (— — 7 2 A) v 2 . (3.48) 
It is now easy to check that 0) satisfies the constraint equation 

[&l + D(A) + XK] 0) = (3.49) 

with cosmological constant A if and only if \I/(z/, 0) satisfies 

[dl + D(A) + K]^(u,^) = (3.50) 

where 

A=y + ^|^(A-1), 0=^0, and = *(»,</>) . (3.51) 

Consequently the two theories are isomorphic. Because of this isomorphism, the gravi- 
tational meaning of the coupling constant A is surprisingly simple: it is related to the 
cosmological constant A. 

Suppose we want to consider the Hamiltonian theory (or the SFM) for zero cosmological 
constant. Then we are interested in the Hamiltonian constraint (3.50) with A = 0. From 
the GFT perspective, on the other hand, the cosmological constant is A which 'runs with 
the coupling constant' A via 3 

A = ^(!- A ) ( 3 - 52 ) 



3 Note incidentally that, contrary to what is often assumed, running of constant under a renormalization 
group flow is not related to the physical time evolution in cosmology [52]. 
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At A = 1, we have A = 0, whence the GFT reproduces the amplitudes of the SFM with 
zero cosmological constant. The question is: What is the space-time interpretation of GFT 
for other values of A? From the perturbation theory perspective, A will start out being zero 
in GFT and, under the renormalization group flow, it will hopefully increase to the desired 
value A = 1. In the weak coupling limit A ~ 0, the SFM will reproduce the amplitudes of the 
theory which has a positive but Planck scale cosmological constant A ps 3/2 , y 2 £ 2 . This is just 
what one would expect from the 'vacuum energy' considerations in quantum field theories 
in Minkowski space-time. As the coupling constant A increases and approaches the SFM 
value A = 1, the cosmological constant A decreases. Now, suppose that the renormalization 
group flow leads us close to but not all the way to A = 1. If we are just slightly away from 
the fixed point A = 1, the cosmological constant A would be small and positive. These 
considerations are only heuristic. But they suggest an avenue by which a fully developed 
GFT could perhaps account for the smallness of the cosmological constant. 



IV. DEPARAMETERIZED FRAMEWORK 

In this section we will use the deparameterized framework which emphasizes the role 
of as internal time. As explained in section II, now we can work in the Schrodinger 
picture, regarding v as the configuration variable and as time. The physical states are 
now represented as functions ^{y) with a finite norm, 

u=4n£ 

and they evolve via Schrodinger equation: 

-idt 0) = VQV{v, 0) = HV(v, 0) . (4.2) 

In contrast to section III, in this section we will not be interested in the kinematical Hilbert 
space or the group averaging procedure. The primary object of interest will rather be the 
transition amplitude 

A(u f , V ;u t ,0) = (u f \e tH ^\u l ) (4.3) 

for the initial physical state \vj) at time 0j = to evolve to \vf) at time 0j = ip. From 
our discussion in section II, one would expect this amplitude to equal the physical scalar 
product ([z//,y?] + , [^j,0] + ) p hy = G(i/f, ip; 0) considered in section III. This is indeed the 
case. For, the positive frequency solution ^ Ui ^ t = [i^,0j] + obtained by group averaging the 
kinematic basis vector |z^,0j) is given by 

0) = fdk (e fe (^) e-*"**) e^W e k (u) (4.4) 

(see Eq.(2.6)) so that the physical scalar product between positive frequency solutions 
[vi,<f>i]+ and [^/,0/]+ is given by 

([^,0/]+, h,0 i ] + ) phy = Jdke^f-^ e k (vi)e k (v f ) (4.5) 

(see Eq (2.9)). The right hand side is precisely the expression of the transition ampli- 
tude (vf\ e lHip \vi) = fdk (uf\ e lHip \k)(k\vi) . Since e k {y) = (u\k), we have the equality: 
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G(uf, tp; i/i, 0) = A(uf, tp; v^, 0). However, the interpretation now emphasizes the physical 
time-evolution in <fi generated by H whence A(i/f, tp; z/j, 0) has the interpretation of a physi- 
cal transition amplitude. Therefore, we can literally follow — not just mimic — the procedure 
Feynman used in non-relativistic quantum mechanics [6]. This will again lead to a vertex 
expansion but one which, if terminated at any finite order, is distinct from that obtained in 
section III. 

In spite of important conceptual differences, the mathematical procedure used in this sec- 
tion is completely analogous to that used in section III. Furthermore, this deparameterized 
framework was discussed in greater detail than the timeless framework in [1]. Therefore, in 
this section we will present only the main steps. 

A. Sum over histories 

Following Feynman, let us divide the time interval (tp, 0) into N equal parts, each of length 
e = tp/N , and express the transition amplitude A(vf, tp; 0) as a sum over discretized paths 
a N = (u f = v N , Pjv-i, . . . , Pi, P = Vi): 

A(is f , tp; u h 0) = ^ A(ctn) with A(a N ) = U PnPn _ 1 U Pn _ iPn _ 2 . . . U P2Pl U PlPo (4.6) 

where now U Pn+lPn = (u n+ i\e lfLH \v n ) . The structure of Eq (4.6) parallels that of Eq (3.9) 
in section III A. However, the mathematical 'time interval' At = 1 in section III A is now 
replaced by the physical time interval (tp, 0) and the mathematical 'Hamiltonian' a9 by 
the physical Hamiltonian H = \^d. Furthermore we no longer split the amplitude into a 
gravitational part and a scalar field part and the group averaging parameter a will never 
appear in this section. 

As in section III A, the next step is to make a convenient rearrangement of this sum, 
emphasizing volume-transitions, rather than what happens at each point 0„ = ne of the 
skeletonized time interval. Thus, we first recognize that the volume could remain constant 
for a number of time steps and consider histories crff with precisely M volume transitions 
(where M < N): 

a N = { ( v m, v m-\i ■ ■ ■ , v \, vq); (N m , N m _i, . . . , N 2 , Ni) }, v m ^ v m -i, 

(4-7) 

where um, ■ ■ ■ ,vq denote the volumes that feature in the history <rj)f and Nk denotes the 
number of time steps after which the volume changes from Vk-i to v^. The probability 
amplitude for such a history aff is given by: 

A(cn ) = [Uv m vm] N Nm 1 Uvmvm-i ■ ■ ■ [Uvivi] N2 Nl 1 U UlUo [U UoUo ] Nl 1 . (4.8) 

As in section III A, we carry out the sum over all these amplitudes in three steps. First 
we keep the ordered set of volumes (vm, . . . ,vq) fixed, but allow the volume transitions to 
occur at any value = ne in the interval X, subject only to the constraint that the m-th 
transition occurs before the (m+l)-th for all m. The sum of amplitudes over this group of 
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histories is given by 

N-l n m -i N 2 -l 
A N (u M ,...,u )= E ••• E A ^n)- (4-9) 

N M =M N M _ 1 =M-1 N 1= l 

Next we sum over all possible intermediate values of v m such that u m 7^ v m -\i keeping 
u = Ui, um = Vfi to obtain the amplitude A(M) associated with the set of all paths in 
which there are precisely M volume transitions: 

A N {M)= Y, A N (v M ,...,v ) (4.10) 

VM-\,—,V\ 

Finally the total amplitude A(uf, 0; z/j, 0) is obtained by summing over all volume transitions 
that are permissible within our initially fixed skeletonization with N time steps: 

N N 

Myf, <r, »i, °) = E a n( m ) = E E Mvm,...,v )]. (4.ii) 

M=0 M=0 vm-i,—,vi 

As in section III A, since A(vf,ip; i^,0) = {vf\e lHif \vi) , the value of the amplitude (4.11) 
does not depend on N at all; the skeletonization was introduced just to express this well- 
defined amplitude as a sum over histories. Thus, while the range of M in the sum and 
the amplitude A N (M) in (4.11) both depend on N, the sum does not. We can get rid of 
the skeletonization altogether by taking the limit as N goes to infinity, to express the total 
transition amplitude as a vertex expansion in the spirit of the timeless framework of spin- 
foams. Reasoning analogous to that in Appendix A shows that the limit does exist. In this 
limit the reference to the skeletonization of the time interval disappears and volume changes 
can now occur at any time in the continuous interval (0j = 0, 0/ = if). The contribution 
Am from paths with precisely M volume changes has a well defined 'continuous time' limit 
and the total amplitude is given by a discrete sum over M: 

00 

A(u f ,(p; u t ,0) = A M^f,V, ^,0) (4.12) 

M=0 

where the partial amplitudes Am are given by 

A M (vf,(p; Vi,0) = E A(vf,v M -i,---Vi,Vi,y) (4.13) 

VM-l,--;Vl 



^ ] Hv m vm-\Hvm-\v M -1 ■ ■ ■ Hu.^H, 



q \ n k -l V 



( n k 1)' \dHw k w k J T[j^m(HwmWm H WjWj ) 



As one might expect, the final expression involves just the matrix elements of the Hamilto- 
man H = Ve. These are calculated in Appendix C. 
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Thus, the total transition amplitude has been expressed as a vertex expansion (4.12) a 
la SFMs. We provided several intermediate steps because, although the left hand sides are 
equal, the final vertex expansions is different from that obtained in section III A: While 
(4.12) features matrix elements of H = \/0, (3.27) features matrix elements of itself. The 
existence of distinct but equivalent vertex expansions is quite surprising. In each case we 
emphasized a distinct aspect of dynamics: the timeless framework and group averaging in 
(3.27), and relational time and deparametrization in (4.12). 

B. Perturbation expansion 

This vertex expansion can also be obtained as a perturbation series that mimics GFTs. 
As in section III, the perturbative approach avoids skeletonization altogether and has the 
advantage that it guarantees a convergent series. Furthermore, since this deparametrization 
approach does not refer to an integral over a, the assumption of interchange of the integral 
and the sum over M that was required in section III B is no longer necessary. 

Let us now focus on the Hamiltonian operator H = \/B (rather than on used in 
section III B) and decompose it into a diagonal part D and the remainder, non-diagonal 
part K which is responsible for a volume change. Finally, let us set H\ = D + XK where A 
will serve as a marker for powers of K, i.e., the number of volume changes in the expansion. 
Then, by working in the appropriate interaction picture, we obtain: 

oo 

A x (is f , <p; Ui, 0) = J2 * M A M (v f , <p; "i, 0) (4.14) 

M=0 

where A M is again given by (4.13). This power series in A is reminiscent of what one finds in 
GFTs. If we set A = 1 at the end of this derivation, we recover the vertex expansion (4.12) 
a la SFMs. For a discussion of the intermediate steps, see [1] and Appendix A. 

C. Satisfaction of the Schrodinger Equation 

Recall that in the deparametrization scheme, the Schrodinger equation (4.2) incorporates 
both the quantum constraint and the positive frequency condition. By its very definition, 
the exact transition amplitude A(i/f, <p; Vi, 0) satisfies this Schrodinger equation. As a check 
on the perturbative expansion (4.14) we are led to ask whether the Schrodinger equation 
would be satisfied in a well-controlled approximate sense if we were to truncate the series 
on the right side of (4.14) at a finite value, say M* of M. We will now show that this is 
indeed the case. 

Since H\ = D + XK, the schrodinger equation would be solved order by order in pertur- 
bation series if for each M we have: 



(id v + D f ) A M (u f , p; Vi, 0) + KfA M _i(uf, tp; v h 0) = . 



(4.15) 
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Using the expression of the partial amplitudes A M we are then led to ask if 



...,V X , Vi\ if) + KfA M -l(Vf, VM-2, v%\ f) 



V M-1, — , V \ Vm-1 



(4.16) 



vanishes for each M. Using the expression (4.13) of A{i/f, vm-i, ■ ■ ■ v±, (p), one can 
readily verify that this is indeed the case. As in section III C, the equation is satisfied 
'path by path', i.e., already by the intermediate amplitudes A(i/f,VM-i,---Vi,Vi', ( p) and 
A{v h v M -i,- ■ - V^Vi^). 

Thus we have shown that the vertex expansion resulting from the perturbation series 
satisfies quantum dynamics in a well-controlled fashion: If we were to terminate the sum at 
M = M*, we would have 



This brings out the precise sense in which a truncation to a finite order of the vertex 
expansion incorporates the quantum dynamics of the deparameterized theory approximately. 



Because LQC is well-developed in the Hamiltonian framework, it provides an interesting 
avenue to probe various aspects of the spin foam paradigm. For definiteness we focused 
on the Friedmann model with a massless scalar field as source. We used the group av- 
eraging procedure that is available for general constrained systems as well as the natural 
deparametrization, with as the emergent time variable, that is often employed in LQC. 

Group averaging provides a Green's function G(uf, 0/; 0j) representing the inner prod- 
uct between physical states extracted from the kinematic kets \vf,(f)f) and \ui,fa). The 
Schrodinger evolution of the deparameterized theory provides the transition amplitude 
A(vf,<f>f;i/i,<f>i) for the physical state \ui) at the initial instant fa to evolve to the state 
\uf) at the final instant of time 0/. We saw in section IV that the two quantities are equal. 
But they emphasize different physics. Following the general procedure invented by Feyn- 
man to pass from a Hamiltonian theory to a sum over histories, we were able to obtain a 
series expansion for each of these quantities — Eq (3.27) for G(vf,4>f, Vi,fa) and Eq (4.12) 
for A(i/f,(j)f; u i: fa) — that mimic the vertex expansion of SFMs. In section III, we had to 
make one assumption in the derivation of the vertex expansion of G(iSf, 0^; i/j, fa): in the 
passage from (3.35) to (3.36) we assumed that the integration over a of the group averaging 
procedure commutes with an infinite sum in (3.35). Since the integration over a is by-passed 
in the deparameterized framework this assumption was not necessary in our derivation of 
the vertex expansion of A(uf, <pf \ z/j, fa) in section IV. 

Detailed parallels between our construction and SFMs are as follows. The analog of the 
manifold M with boundaries Si, Sf in SFMs is the manifold Vxl, where V is the elementary 
cell in LQC and X, a closed interval in the real line (corresponding to r e [0, 1] in the timeless 
framework and G [<f>f,fa] in the deparameterized). The analog of a triangulation in spin- 
foams is just a division of V x I into M parts by introducing M — 1 time slices. Just 



M* 



(id v + D f + \K) [ ^ M A M (v f , v h 0)] = <D(\ M * +1 ) 



(4.17) 



M=0 



V. DISCUSSION 
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as the triangulation in SFMs is determined by the number of 4-simplices, what matters in 
LQC is the number M; the precise location of slices is irrelevant. The analog of the dual- 
triangulation in SFMs is just a 'vertical' line in V x X with M marked points or 'vertices' 
(not including the two end-points of X). Again, what matters is the number M; the precise 
location of vertices is irrelevant. Coloring of the dual-triangulation in SFMs corresponds to 
an ordered assignment (vm, vm-i, ■ ■ ■ vi, ^o) of volumes to edges bounded by these marked 
points (subject only to the constraints vm — Vf, v = z/j and u m ^ v m -i). Each vertex signals 
a change in the physical volume along the quantum history. 4 The probability amplitude 
associated with the given coloring is given by A(uf, . . . , u ; <pf, (pi) in the group averaging 
procedure (see Eq (3.26)) and by A(uf, . . . , u ; <p) in the deparametrization procedure (see Eq 
(4.13)). A sum over colorings yields the partial amplitude associated with the triangulation 
with M 'vertices'. The Green's function G(uf,(j)f] z/j,0j) and the total transition amplitude 
A(vf, if] z/j, 0) are given by a sum over these M-vertex amplitudes. 

Thus, the physical inner product of the timeless framework and the transition ampli- 
tude in the deparameterized framework can each be expressed as a discrete sum without 
the need of a 'continuum limit': A countable number of vertices suffices; the number of 
volume transitions does not have to become continuously infinite. This result supports the 
view that LQG and SFMs are not quite analogous to quantum field theories on classical 
space-times. Discrete quantum geometry at the Planck scale makes a key difference. In 
sections III B and IV B we were able to obtain the same vertex sum using a perturbative 
expansion, in a coupling constant A, that is reminiscent of GFTs. In sections III C and 
IV C we showed that this is a useful expansion in the sense that the Green's function and 
the transition amplitude satisfy the dynamical equations order by order in A. Thus, if we 
were to truncate the expansion to order M, the truncated Green's function and transition 
amplitude would satisfy the dynamical equations up to terms of the order 0(A M+1 ). Finally 
in section III D we showed that the coupling constant A inspired by GFTs is closely related 
to the cosmo logical constant. This interpretation opens a possibility that a detailed study of 
the renormalization group flow in GFT may be able to account for the very small, positive 
value of the cosmological constant. 

Taken together, these results provide considerable concrete support for the general 
paradigms that underlie SFM and GFT. 5 However, we emphasize that this analysis has 
a key limitation: We did not begin with a SFM in full general relativity and then arrive 
at the LQC model through a systematic symmetry reduction of the full vertex expansion. 
Rather, we began with an already symmetry reduced model and recast the results in the 
spin foam language. Reciprocally, a key strength of these results is that we did not have to 
start by postulating that the physical inner product or the transition amplitude is given by 
a formal path integral. Rather, a rigorously developed Hamiltonian theory guaranteed that 
these quantities are well-defined. We simply recast their expressions as vertex expansions. 
The complementarity of the two methods is brought to forefront in the recent work [54] on 
spin-foams in the cosmological context. There, one begins with general spin foams, intro- 



4 In the Bianchi models there are additional labels corresponding to anisotropics [53] . These are associated 
with the faces of the dual graph, and are thus analogs of the spin labels j associated with faces of general 
spin foams. 

5 But it also brings out the fact that the term 'third quantization' that is sometimes used in GFTs can be 
misleading in other contexts. In cosmology, the term is often used to signify a Fock space of universes, 
where the 'single universe sector' is described by the theory presented here. 
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duces homogeneity and isotropy only as a restriction on the boundary state and calculates 
just the leading order terms in the vertex expansion. By contrast, in this work we restricted 
ourselves to homogeneity and isotropy at the outset but calculated the physical inner prod- 
uct (or, in the deparameterized picture, the transition amplitude) to all orders in the vertex 
expansion. 

It is often the case that exactly soluble models not only provide support for or against 
general paradigms but they can also uncover new issues whose significance had not been 
realized before. The LQC analysis has brought to forefront three such issues. 

First, it has revealed the advantage of adding matter fields. It is widely appreciated 
that on physical grounds it is important to extend SFMs beyond vacuum general relativ- 
ity. However what was not realized before is that, rather than complicating the analysis, 
this generalization can in fact lead to interesting and significant technical simplifications. 
This point is brought out vividly by a recent analysis of Rovelli and Vidotto [55]. They 
considered a simple model on a finite dimensional Hilbert space where there is no analog 
of the scalar field or the possibility of deparametrization. There, individual terms in the 
vertex expansion turn out to be well defined only after a (natural) regularization. In our 
example, the presence of the scalar field simplified the analysis (in the transition from (3.24) 
to (3.26)) and individual terms in the vertex expansion are finite without the need of any 
regularization. Furthermore, this simplification is not an artefact of our restriction to the 
simplest cosmological model. For example, in the Bianchi I model the Hamiltonian theory 
is also well-developed in the vacuum case [56]. Work in progress by Campiglia, Henderson, 
Nelson and Wilson-Ewing shows that technical problems illustrated in [55] arise also in this 
case, making it necessary to introduce a regularization. These problems simply disappear 
if one also includes a scalar field. A qualitative argument suggests that the situation would 
be similar beyond cosmological models as well. 

Second, it came as a surprise that there are two distinct vertex expansions: Group 
averaging provides one that mainly uses the matrix elements of 6 while the deparameterized 
framework provides one that uses only the matrix elements of y/Q. This is not an artefact 
of using the simplest cosmological model. Work in progress indicates that the situation is 
similar in the anisotropic Bianchi models. Indeed, from a Hamiltonian perspective, it would 
appear that distinct vertex expansions can arise whenever a well-defined deparametrization 
is available. This raises an interesting and more general possibility. Can there exist distinct 
spin foam models — constructed by using, say, distinct vertex amplitudes — for which the 
complete vertex expansions yield the same answer? Finite truncations of these expansions 
could be inequivalent, but each could be tailored to provide an excellent approximation to 
the full answer for a specific physical question. One may then be able to choose which 
truncated expansion to use to probe a specific physical effect. 

The third issue concerns three related questions in the spin foam literature: i) Should the 
physical inner products between states associated with spin networks be real rather than 
complex [31]? ii) In the classical limit, should one recover cos S in place of the usual term 
e lS , where S is the Einstein Hilbert action [32, 33]? iii) Should the choice of orientation 
play a role in the sum of histories [49]? In the LQC example we studied in this paper, these 
three questions are intimately related. The inner product between the physical states [is, <p] + 
determined by the kinematic basis vectors — which are the analogs of spin networks in this 
example — are in general complex (see Eq (3.27)). However, if we had dropped the positive 
frequency requirement, the group averaged inner products would have been real (see Eq 
(3.24)). The situation with action is analogous. And, as we show in the next paragraph, 
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the positive frequency condition also selects a time-orientation. 

Since this is an important issue, we will discuss it in some detail. Let us begin with 
the classical theory. The phase space is 4-dimensional and there is a single constraint: 
C(v, b; (f),p^) := Gp^Sn (ip^ 2 ) b 2 = 0. Dynamics has two conceptually interesting features. 
First, given a solution (u(t),(p(t)) to the constraint and dynamical equations, (— z/(t), 4>{t)) 
is also a solution (where t denotes proper time). They define the same space-time metric 
and scalar field; only the parity of the spatial triad is reversed. Therefore (z/(t),0(t)) — >■ 
(— v(t), 4>(t)) is regarded as a gauge transformation. The second feature arises from the 
fact that the constraint surface has two 'branches', p^ > and p^ < 0, joined at points 
p^ = which represent Minkowski space-time. As is usual in quantum cosmology, let us 
ignore the trivial flat solution. Then each of the two portions of the constraint surface 
defined by the sign of p^ is left invariant by dynamics. Furthermore, there is a symmetry: 
Given a dynamical trajectory (z/(t),0(t)) in f + , there is a trajectory (v(t), —<f>(t)) which 
lies in T~ . This represents a redundancy in the description in the sense that we recover all 
physical space-time geometries g a b{t) even if we restrict only to one of the two branches T ± . 
In particular, the dynamical trajectories on r + , for example, include solutions which start 
with a big-bang and expand out to infinity as well as those which start out with infinite 
volume and end their lives in a big crunch. The difference is in only in time orientation: If 
we regard <fi as an internal or relational time variable and reconstruct space-time geometries 
from phase space trajectories, space-times obtained from a trajectory on f + defines the 
same geometry as the one obtained from the corresponding trajectory on f _ but with 
opposite time orientation. As in the Klein-Gordon theory of a free relativistic particle, 
this redundancy is removed by restricting oneself either to the p^ > sector or to the 

< sector. In the quantum theory, then, the physical Hilbert space is given by solutions 
ty(v,(j)) to the quantum constraint (2.2) which in addition have only positive (or negative) 
frequency so that the operator p^ is positive (or negative) definite. (They are also invariant 
under parity, \I/(z/, 0) = v, </>)). Thus, the LQC example suggests that in general SFMs 
one should fix the time-orientation, lending independent support to the new ideas proposed 
in [49]. Reality of the physical inner products between spin network states [31] and the 
emergence of cos S in place of e lS [32, 33] can be traced back to the fact that in most of the 
SFM literature one sums over both orientations. However, our analysis provides only a hint 
rather than an iron-clad argument because all our discussion is tied to LQC models where 
symmetry reduction occurs prior to quantization. 

We conclude with an observation. We have recast LQC as a sum over histories. However, 
this is different from a Feynman path integral in which the integrand is expressed as e tS , 
for a suitable action S. This step was not necessary for the goals of this paper. However, it 
is of considerable interest, especially in the cosmological context, for certain physical issues 
such as the emergence of the classical universe and semi-classical corrections to the classical 
theory. Such a path integral formulation of LQC does exist [57] and will be discussed 
elsewhere. 6 



A path integral formulation of polymer quantum mechanics was carried out independently by Husain and 
Winkler [58]. 
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Appendix A: Limit in Eq (3.18) 



It is convenient to rewrite An(um, ■ ■ ■ , ot) defined in (3.13) in the following way: 



r] Wu M VM ■ 
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Uv M -l"M-l 
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N-l N M -1 N 2 -l Vtt -|AT m Vtt -| JVi 

E E -E 

N M =M N M _ 1= M-1 Ni=l 



A(u M , ...,v ;a) = (-ia) M Qu M u M ^ 



6 e 



(Al) 



Our aim is to calculate the limit N — > oo of (Al) and show that is given by A(vm, ■ ■ ■ , W, 
of Eq (3.18) which we rewrite as 



/;dr M /?'dr M _ 1 ...J7dr 1 e^...e^ 

where 

We start by calculating the N 3> 1 behavior of the terms appearing in (Al). These are: 

U Wm = -^Qu m+1 u m + 0(N~ 2 ), (A4) 



(A2) 
(A3) 



\U liV = N\o g U l , MVM 

\yvMVM\ c 

= e N (~^M^M+o(N- 2 )) 

= e -i<x6» M »M + o(N~ 1 ), (A5) 
[U VMVMt ...U vavtt ]~ 1 = 1 + OiN- 1 ), (A6) 



_ e A r m(logl/„ m _ 1I , m _ 1 -logt/ I/m „ m ) 

= e N m (b m /N+0(N- 2 )) 

= e^ bm +0(N m N- 2 ), (A7) 

with b m given in (A3). In (A5) and (A7) we have used the fact that the multivalued nature 
of the log function does not affect the final result: e N ( lo & x + 2mk ) — e Niogx w ] iere ^ e ^ reflects 
the multiple values that log can take. 
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We now substitute expressions (A4) to (A7) in (Al) to obtain 

A N (u M , ...,u ;a) = [{-ia) M Q VMVM _, . . . Q VlV0 e~^^ M N ~m + o^"^ 1 )] x 



M 

n 

m=l 



N m+ i-l 



J2 e^ b - + 0{N m N' 2 ) 



N m =m 



(A8) 



where the product denotes the M nested sums in (Al). Each sum in (A8) has two 



terms. The first one gives a contribution of ^2 Nm e" 



O(N) while the second one is 



X)jv m 0(N m N- 2 ) ~ 0(1). The M sums then give a contribution of order [0(N) + 0(1)] M ~ 
0(N M ) + 0(N M ~ 1 ). By combining this with the first factor of (A8), we find that the 
non-vanishing contribution comes from the first terms of the sums: 



A N (v M ,...,v ;a) =(-ia) 
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N m =m 



e N 



+ 0(N- 1 ). 



(A9) 



Eq (A9) has all the pre-factors appearing in (A2). It then remains to show that N~ M 
times the sums in (A9) limits to the integrals in (A2). But this is rather obvious, as the 
sums can be seen as Riemann sums for the integrals. Specifically, 



lim N~ M 
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(A10) 



where, in the second line, we have slightly changed the limits on the sums, introducing an 
0(iV _1 )-term which vanishes in the limit. This concludes the proof of the limit (3.18). 



Appendix B: General Integrals in Eq (3.18) 

The integrals over r appearing in the amplitude for a single discrete path (3.18) can be 
evaluated for a general sequence of volumes (u M , v ) with the result given by (3.21). In 
this appendix we will perform these integrals first for the case where all are distinct and 
then for the general case. The amplitude for a single discrete path given by (3.18) and (3.19) 
is 

A(u M , . . . , i/o, a) =J^dr M J^dm^ . . . J^ drie ^(Ar-r M ) a e^ M (-i a e UMVMi ) x 

e -K^M-TM-i)aB VM _ lVM _ 1 ^ e -i(r 2 -Ti)ae nn (—ja© ) e i7lae D1i (Bl) 
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This expression can be written in terms of the following integral. 

I(x M , ...,x , At) = J^ T dr Af J T M dr M -i ■ ■ ■ Jo^i(i) M e^-™*" e ^M-T M -,)x M -, ( B2 ) 



We will first evaluate this integral for the case where all Xi are distinct. By induction on 
M — the number of vertices or the number of times that x changes value — we will show 
that when the Xj are all distinct the integral is given by 



M 



I(x M ,...,x , At) = J2 u m , — 7 ( B3 ) 

i=0 lljfr\ X i X 3l 

This is true by inspection for M — 0. If we assume that (B3) holds for M we can evaluate 
the integral with M + 1 vertices. 

I(x M+ i,x M , ...,x , At) = J^ r dT M+1 ie^ AT -™+^ x ^I(x M , ...,x , t m +i) (B4) 



Y\j 1 ii{ x i~ x 3) 



M p ix,Ar M 



= 2. 



j^o n?^ ( x « i=o rif^i ^j) 

In the first step we recognized that the M + 1-th integral contains the M-th and then, in 
the second step, we inserted the assumed result for the M — th integral. In the second step 
the integral over t m+1 is carried out. Finally using the identity 

M+l ^ 

E ftm+T/ _ ; = ( B5 ) 

i=l LLjfr \ X 'i X 3> 



The integral can be written as 



M + 1 e ix t Ar 



I{x M+ i, x M ,..., x , At) = ^ i-rM+i/ — - — 7 ( B6 ) 

i=0 iij^i \ X i x i) 

Therefore if (B3) holds for M it also holds for M + l, thus by induction it holds for all 
M > 0. 

If the Xi are not distinct, if there exist i,j such that Xj = Xj, then the proof follows 
in a similar way. The key element is that the integral /(%,...,xo) is independent of the 
order of the x^s. This can be seen by rewriting the integral in terms of the time intervals 
ATj = Tj+i — Tj where t = and r m+ i = At. 

At) = /^ r dAT M dAT M -i..dATo 5(AT m + ... + At - At) (B7) 

(l) M e i A ™ XM e i A ™~ 1 XM ~ 1 e * An x i e i Ar xo 

It is clear that this is symmetric under the interchange of Xi with Xj for all i,j, so the integral 
is independent of the order of the sequence Xi. Since the integral is independent of the order 
of the values X; t it should be characterized by the distinct values, labeled by and their 
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multiplicity rij. Where n\ + . . . + n p — M + 1. Given a set of values Xi we will evaluate the 
integral for the case where they are organized such that any x^ sharing the same value are 
grouped together. Doing so the integral simplifies to 

I(y p ,n p ,...,y 1 ,n 1 ,Ar) = /^dr M /™dr M _i .. ./JdnWV^^-i)* (B8) 

e i(r ni +...+n p _ 1 -Tn 1 +...+n p _ 2 )y P -l g*(T ni +„ 2 -T ni )V2 ^T nx J/1 



By induction on p, the number of distinct values, we show that this integral is given by 



I(y p ,n p , ...,y 1 , ni, At) 



9 



(np - 1)! V% 



fn fc - 1)! 



n v — l 



1 



9 



ni-1 P 



JyiAr 



■■■(m-l)! UJ 

(B9) 
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n fc -l P 



E FT'' 



(n fc - 1)! V% 
For p — 1 (B8) can be easily evaluated giving 

/(y 1 ,n 1 ) = /^dr Bl _ 1 ...J7dr 1 (i) B1 - 1 e i « 



-IJyxAr _ (»Ar)"l 1 tmAr 



(ni-1)! 



(BIO) 



Til — 1 



3 »J/»At 



(ni - 1)! 



If we assume that (B9) holds for p distinct values then we can evaluate it for p + 1 distinct 
values as follows. 



I(y p+ i, n p+1 , y p ,n p ..., y u m, At) = J^ T dr M . . . jT np+L+2 dT M -n p+1 +i 

W n P+1 -l e i(Ar-r M -„ p+1+I fe +1/(2/p) ^ _ _ _ j ^ ^ rM _ np+i+1 ) 



(Bll) 



Plugging in the assumed result for p distinct values and performing the integrals over r we 
obtain 



HVp+i, n p+1 , ...,y 1 ,n 1 ,Ar) = Yl , _ ^, 



L (n fc - 1)! \9j/ fc 



. , \ n k -l P 

" 1 £ 



g e*+ lAr (iAr)™^ 1 - 



(B12) 



(j/i - y P +i) np+1 ^ (j/i - y P +i) m K+i - m)! 

We recognize that the term in brackets can be written as derivatives with respect to y p+ i of 
a simple function. 



V 1 / q \n k -l P 

I(y p+1 , n p+1 , y p ,n p ..., y u m, At) = JJ ^— _ ^ J ^ 



(n p+ i - 1)! V^P+i 



9 



(B13) 
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Finally simplifying the expression and using eqn (B5) we obtain 

p +1 y ( d \ ™ fe_1 P+1 e iy%! ^ T 

-,«.,».. at) =n 5^ (g^j gn^fe^) < B14 > 

Thus if (B9) holds for p then it also holds for p + 1, so it is true for all p > 0. Using this 
result we find that the contribution due to each discrete path is 

A(v M ,...,v ,a) = (O )(0 )...(0^ 1 )(0, 1 ,o) (B15) 

11 Ufl ) 



l (n fc - 1)! v 96 ^/ lTy©^ - 



where Wi label the distinct values taken by v along the path and rii the multiplicity of each 
value. 



Appendix C: Eigenstates and Operator functions of 

In the timeless framework of section III, the vertex expansion mostly featured matrix 
elements ® VmVn = {v m \®Wn)- These are easy to evaluate directly from the definition (2.3) of 
0. In the deparameterized framework of section IV, on the other hand, the vertex expansion 
involves matrix elements of \^Q. To evaluate these one needs the spectral decomposition of 
0. In the first part of this Appendix we construct eigenstates of and discuss their relevant 
properties. In the second part we use these eigenstates to evaluate the matrix elements 
functions of 0, including v^©- 



1. Eigenstates of 

Recall that is a positive, self-adjoint operator on "Hf[„ v . By its definition (2.3), it follows 
that preserves each of the three sub-spaces in the decomposition "Hf[„ v = H- © "Ho © 
spanned by functions with support on v < 0, v = and v > respectively. In particular, 
\v — 0) is the unique eigenvector of 0, with eigenvalue 0; H is 1-dimensional. Our first 
task is to solve the eigenvalue equation for a general eigenvalue uj\: 

ee k (u)=co 2 k e k (u). (CI) 

This task becomes simpler in the representation in which states are functions x(b) of the 
variable b conjugate to v. 7 

m:= Jk^ e i*m. (C2) 

V 71 AC 

v=Anl 



7 Our normalization is different from that in [39]. The wave function ^(v) in [39] is related to the one here 
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In this representation, the eigenvalue equation (CI) takes the form of a simple differential 
equation 

(6x*) (6) = -127TG (^d b y X k(b) = ul Xk {bl (C3) 



whose solutions are 

X k(b) = A(k) e tklog{tan ^ with u 2 k = 12txG k 2 , (C4) 



where A(k) is a normalization factor and k G (—00, 00). k = yields a discrete eigenvalue 
Uk = and in the z/ representation the eigenvector can be expressed simply as e (u) = 8 0;l/ . 
Eigenvectors with non-zero eigenvalues can also be expressed in the v representation by 
applying the inverse transformation of (C2) to (C4): 



e k {y) = A{k)^^fl /e °db e -^b e ikio S (t a n^) where k ^ o _ 



(C5) 



Let us note two properties of these eigenvectors. First, and e__fc have the same eigen- 
value and so the wj|;-eigenspace is two-dimensional. Second, the vectors e^{y) we have ob- 
tained have support on both v > and v < 0. However, since O preserves the sub-spaces 
H±, it is natural to seek linear combinations e^(z/) of e^iy) and e_fc(i/) which lie in these 
sub-spaces. In particular, this will simplify the problem of normalization of eigenfunctions. 

Let us begin by rewriting the integral in (C5) as a contour integral in the complex plane. 
Recalling that v = Ai Q n and setting z = e Mo we obtain 



where C is the unit semicircle in counterclockwise direction in the upper half, Qz > 0, of the 
complex plane. As remarked earlier, e^iy) = A(k)^ l \v\/ir J(k, v/M ) has support on both 
positive and negative values of v — 4£ n. Now, the second independent eigenfunction e-k{y) 
with the same eigenvalue uj\ can be represented in a similar fashion by setting z = —e Mo . 
The result is a contour-integral along the unit semicircle in counterclockwise direction in 
the lower half, Qz < of the complex plane. By combining the two integrals, we obtain a 
closed integral along the unit circle: 



2 ^ j z- 2 ^ 1 (Lpf) ^ dz = l - {e^ 2 J(k, n) + e~^ 2 J(-k, n)) =: I(k, n) . (C7) 



Being a linear combination of e^iy) and e^[y\ this I(k,n) gives also an eigenfunction of 
with eigenvalue io\. Moreover, using elementary complex analysis, one finds that it has 
support only on positive n: 



1 d 2n 



I(k,n) = < ( 2 ") ! ^ s =o Kl+sJ (C8) 
n < 0. 
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Repeating the argument but taking z = e Mo and z = — e Mo one obtains 

\ (e-^ 2 J(k, n) + e^ 2 J(-k, n)) = / z 2n ~ x dz = I(k, -n) 

which has support only on negative n. Thus, the basis we are looking for is given by 



(C9) 



eJW == \ (e ±7Tk/2 e k (u) + e^ 2 e. k (u)) = A{k) ^ 



(CIO) 



By construction, G "H±. 

Next, let us calculate the normalization of these vectors. It is convenient to introduce 
kets \k±) such that (y\k±) = e^(u). Then, it is clear that (k' ± \k=f) = 0. To calculate 
the nontrivial inner product, (k' ± \k±), let us return to the b representation. There, the 
functions describing the states \k±) are 



Xt(b) = 



A(k) 



fe/2 e ifelog(tan%5) eT7rfe /2 e -ifclog(tan%^) 



and their inner product is given by [39] 

(k' ± |*±) = ff°db \A(k)\ 2 ^ k ,(b) \m\ x^b) 



(Cll) 



(C12) 



where \2idf,\ is the absolute value of the volume operator v = 2idb. Simplification occurs 
because e k (u) have support only on positive/negative v values. Because of this property, 
one can replace \d\\ in (C12) by ±<9{,. The calculation now reduces to a straightforward 
integration. The result is 



(k'±\k±) = \A(k)\ 2 27rk sinh(vrA;) 5(k',k). 



(C13) 



2. Matrix Elements for /(9) 

We will now use the eigenbasis | ± k) of to calculate the matrix elements 
{4n£ \f(Q)\4m£ ), of the operators of the form /(0), for a measurable function /. Through- 
out this section, the normalization factor A(k) is chosen to be unity. From the normalization 
condition (C13) with A(k) = 1, we have the following decomposition of the identity: 



Joe 




dfc 



2irk sinh(7rfe) 



(\k+){k+\ + \k-){k- 



(C14) 



which can be inserted in (4n£ \f(Q)\4m£ ). If m and n have different signs, the result is 
zero. It suffices to consider the case where both are positive. By writing (An£ \k+) in terms 
of derivatives (see equations (CIO) and (C8)), one obtains 



(4n£ \f(e)\4m£ ) = 



2 X /Wm d 2m d 2n 



(2n)\(2m)\ ds 2m dt 2n 



7(6) 



s=t=0 



1+sl-t 
l-sl+t 



(C15) 
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with -F/(e) the 'generating function' given by 8 

F m {x) = j^ki^f. (C16) 

We now give the generating function for v^©- It is also useful (at least to check normalization 
factors) to write down the generating functions for operators whose matrix elements are 
known, namely 9 and the identity /. These generating functions are given by, 

F^x) = -2^1og(l+a;) + logr(l/2 + 2^)^ (C17) 

F^(x) = Vl2^G {^f- - 1^(1/2 + i^)) (CIS) 

Fe(x) = 12^(^-^(1/2 + ^)) (C19) 

where T(z) is the gamma function, and ip(z) = T'(z)/T(z) the polygamma function. 

In obtaining these functions, it is useful to observe the following relations among them: 

F^(x) = -iVvhrG x^-F^x) (C20) 

F e (x) = -iVUnGx-j^F^x), (C21) 

which can be derived from (C16). 

We will conclude by noting that the matrix elements for the evolution operator U (tp) = 
e i>fiVe are eaS y t find; From (C16) one sees that Fu^(x) = F/(e^ 127rG</, x). 
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